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Abstract 

We obtain the recursive identities for the joint moments of the traces of the powers of 
the resolvent for Gaussian ensembles of random matrices at the soft and hard edges of the 
spectrum. We also discuss the possible ways to extend these results to the non-Gaussian 
case. 

1 Introduction 

Consider the Gaussian Orthogonal Ensemble (GOE) of real symmetric n x n random matrices 

A n = —= {a,ij) ij=1 , (1) 



where {a^- = ciji}i<j are independent iV(0, 1 + Sy) random variables. GOE is the archetypal 
example of a Wigner real symmetric random matrix where the matrix entries {ay = dji}i<j are 
assumed to be independent up from the diagonal, centralized, and to have the same variance 
(except, possibly, on the diagonal). It follows from the classical Wigner semi-circle law (([20 , 
[21j . pQ) that the empirical distribution function of the eigenvalues of A n converges as n — > oo 
to the limiting distribution with the probability density ^-y/i — x 2 supported on the interval 
[—2,2]. Celebrated work by Tracy and Widom (see [17J for the GOE case) proved that the 
largest eigenvalues of A n deviate from the right edge of the spectrum on the order of n~ 2 / 3 . In 
particular, Tracy and Widom calculated the limiting distribution of the largest eigenvalue: 

JirnP (\ max < 2 + in" 2 / 3 ) = Fi(x) = exp ^-1/2 j q{t) + {t - x)q 2 {t)dt^j , (2) 



*Department of Mathematics, University of California at Davis, One Shields Ave., Davis, CA 95616, USA. 
E-mail address: soshniko@math.ucdavis.edu. Research was supported in part by the NSF grant DMS-0707145. 



1 



where q(x) is the solution of the Painleve II differential equation q"(x) = xq(x) + 2q 3 (x) with 
the asymptotics at infinity q{x) ~ Ai{x) as x —* +00. Here ^4z(x) denotes the Airy function. 

To consider the joint distribution of the largest eigenvalues at the edge of the spectrum, we 
rescale the eigenvalues as 

A W = 2 + e] n) n- 2 / 3 , , 1.2 n. (3) 

where A< n) > \^ . . . > are the ordered eigenvalues of A n . It then follows from the results 
of [13], [18] that the random point configuration {C- , 1 < j ' < n} converges in distribution 
on the cylinder sets to the random point process on the real line with the Appoint correlation 
functions given by 

/ \i/2 
p k (x 1 ,...,x k ) = [detiKix^Xj))^^ , (4) 

where K(x, y) is a 2 x 2 matrix-valued kernel with the entries 

Kn(x,y) = K Airy (x,y) + ^Ai(x) ^1 - J Ai(z)dz^ , (5) 



1 

2" 

+00 / fX /"+00 



Ku(x,y) = -d y K Airy (x,y) - -Ai(x)Ai(y), (6) 



K 21 (x,y) = -J K Airy {z,y)dz + -\^J Ai{z)dz + J Ai{z)dz J Ai(z)dzj,(7) 
K 22 (x,y)=K n (y,x), (8) 
and the Airy kernel K A i ry {x,y) is defined as 

K Ai ry{x,y = / + z)At(y + z)dz = — ^ ^ — — — . 9 

Jo x-y 

Therefore, the fe-point correlation function of the limiting random point process is given by the 
square root of the determinant of the 2k x 2k matrix defined in (j4j) and ([5]). One can also rewrite 
the fc-point correlation function in the pfaffian form (see e.g. [9]) which shows that the limiting 
random point process belongs to the family of the pfaffian random point processes (see e.g. 
|14j). In particular, the right-most particle of this pfaffian random point process is given by the 
Tracy- Widom distribution ([2]). Moreover, it was shown in [12] that the asymptotic behavior of 
the largest eigenvalues is universal for Wigner real symmetric matrices with sub-Gaussian and 
symmetrically distributed entries. 
Define 

G n (z) = (A n -2- zn- 2 / 3 )- 1 (10) 

for complex z with non-zero imaginary part Qz 7^ 0. Here and throughout the paper, we will 
use (A n — z)^ 1 as the shorthand notation for the resolvent matrix (A n — zld)^ 1 . 



2 



Let 



n 

g nAz ) = n-^TrG k n (z) = n-^Tr(A n - 2 - zn 2 l 3 y k = - Z y k (11) 

1 

for positive integers k = 1, 2, . . . . It can be shown that for k > 2, g n ^( z ) is a "local" statistic of 
the largest eigenvalues in the GOE. Indeed, only eigenvalues from the 0(ra~ 2 / 3 )-neighborhood 
of the right edge of the spectrum give non- vanishing contribution to g nj k(z) in the limit n — > oo. 
For example, the joint contribution of the eigenvalues from (— oo, 2 — 6} can be trivially bounded 
in absolute value by n l ~ 2k ^\5+ zn~ 2 / 3 \ k = o(l) for large n uniformly in z with bounded from 
below. More delicate estimates involving the asymptotics of the one-point correlation function, 
imply that the joint contribution of the eigenvalues from (— oo,2 — ra~ 2 / 3+e ) to g n ,k( z ) is still 
negligible for all e > and k > 1. Moreover, the one-point correlation function pi(x) of the 
limiting pfaffian random point process defined in ([3])-([9]) decays super-exponentially at +oo and 
grows proportionally to |x| 1//2 at — oo. Consequently, if £ = j 6 Z} is a random point 

configuration of the limiting pfaffian random process then 

/ + 0O 
\x - z\- k Pl (x)dx < oo (12) 
j 

for any integer k > 2. The integral at the r.h.s. of (112f) diverges for k = 1 which emphasizes the 
fact that g nt i(z) is not a "local statistic" as the main contribution to g n ,i( z ) = n~ 2 / 3 x 
Tr(A n - 2 - zn- 2 / ?J )- 1 comes from the eigenvalues in the bulk. Moreover, it could be shown 
from the asymptotics of the GOE one-point correlation function that 

E (Tr(A n - 2 - zn^ 3 )" 1 ) =-n + 0(n 2/3 ). (13) 

Eventhough the eigenvalues from the bulk of the spectrum give the main contribution to the 
mathematical expectation of Tr(A n — 2 — zn" 2 / 3 )" 1 , their joint contribution to the fluctuations 
of Tr(A n — 2 — zn^/ 3 )^ 1 around its mean is much smaller (namely, it can be shown to be of 
order of constant if one smoothes their contribution by a test function with the support inside 
[-2 + 5,2 

delta].) On the other hand, the largest eigenvalues of A n give smaller (namely, of the order of 
0(n 2 ' 3 )) contribution to the mean of Tr(A n -2-zn~ 2 / 3 y 1 , but they give the main contribution 
to the fluctuations of Tr(A n — 2 — zn" 2 / 3 )™ 1 around its mean. This suggests to consider 

<£,l0) = n~ 2/3 (n + TrG n (z)) = rT 2 ! 3 (n + Tr{A n - 2 - zn" 2 / 3 )" 1 ) (14) 

which is a "local" statistic in a sense that the main contribution to g^ ±(z) comes from the largest 
eigenvalues (i.e. the eigenvalues that deviate from the right edge of the spectrum on the order 
of 0(n- 2 / 3 ).) 
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In Theorem 1.1, we obtain the recursive relations on the joint moments of the local linear 
statistics g^ and g n ,k( z )i k > 2. Let K be a multi-index, K = (k\, . . . , kj), j > 1, with 
the components ki, 1 < I < j, nonnegative integers. The number of components j is not fixed. 
We will denote by mx the corresponding joint moment of i(z) and g n ,k{ z )i k > 2, namely: 

^=iE^(^iW) fci n(^w) fc ^- (15) 

Let e/ denote the multi-index with the l-th component equal to 1 and the other components 
equal to zero. 

Theorem 1.1. Let K be a non-zero multi-index, then the following equation holds: 

m K (z + 0(n- 2 / 3 )) - m^ +2ei (l + 0{n~ 2 ' 3 )) - m K+e2 (l + 0{n- 2 ' 3 )) 

- 2Y,lkim K . ei+ei+2 {l + 0{n~ 2 ^)) = 0{n- 1 ^)m K+ei . (16) 
i>i 

Also, the following "boundary" condition holds: 

z + 0(n- 2 / 3 ) - m 2ei (l + 0(n~ 2 / 3 )) - m e2 (l + 0( n - 2 / 3 )) = 0(n^ 1 / 3 )m ei . (17) 



Remark. We will always assume in (|16p that kimK-ei+ei +2 = if ^ = 0. 
Theorem 1.1 will be proved in the next section. Let us now consider the Gaussian Unitary 
Ensemble (GUE) of Hermitian n x n random matrices. 

A n = -L(a jk )l k=1 , (18) 



n 



where {"Siajk = yiakj}j<k arid {Qajk = — 9afcj}j<fe are i.i.d. iV(0, 1/2) random variables, and 
{au}i<i<n are i.i.d. iV(0, 1) random variables. 

The global distribution of the eigenvalues of A n still satisfies the Wigner semi-circle law in 
the limit n — > oo. The limiting local distribution of the largest eigenvalues of A n was calculated 
by Tracy and Widom in [16]. In particular, 

(t-x)q 2 (t)dt ) , (19) 



Jim P [\ rnax < 2 + xn" 2 / 3 ) = F 2 (x) = exp (- J 



where, as before, q(x) is the solution of the Painleve II differential equation with the same 
asymptotics at infinity. 

Consider the same rescaling at the right edge of the spectrum as in the GOE case, namely 

Af ) = 2 + ef ) n- 2 / 3 , j = l,2,..., (20) 
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where a[" } > \f ] . . . > A& are the ordered eigenvalues of A n . It then follows from the results 
of [16] that the random point configuration {^j ? 1 < i < ^} converges in distribution on the 
cylinder sets to the random point process on the real line with the A;— point correlation functions 
given by 

p k (xi,...,x k ) = det(K(x i ,x j )) 1 ^ ij ^ k (21) 

where K(x,y) = Kjn ry (x, y) is the Airy kernel defined in ([9]). The limiting random point process 
belongs to the class of determinantal random point processes (see [IS], [6]). 

Let us use the same notations G n (z), g n ,k(z), g^ ll (z), and Mjc in the GUE case as they 
were defined in (TlO|) . (fTTI) . (fT4"j) . and (fl5|) in the GOE case above. The following analogue of the 
Theorem 1.1 holds: 

Theorem 1.2. Let K be a non-zero multi-index, then the following equation holds: 

m K (z + 0(n" 2 / 3 )) - m K+2ei (l + 0(n- 2 / 3 )) 

- lkim K -e l+ei+2 (l + 0(n" 2 / 3 )) = 0(n-y 3 )m K+ei . (22) 
i>i 

Also, the following "boundary" condition holds: 

z + 0(n- 2 / 3 ) - m 2ei (1 + 0(n- 2 / 3 )) = 0(n- 1 / 3 )m ei . (23) 

We now turn our attention to Wishart (a.k.a Laguerre) ensembles of random matrices. Again, 
we start with the real case. Let A = A n ^ = A=(aij) be a rectangular n x N matrix with 
{ciy, 1 < i < n, 1 < j < N, } real i.i.d. A^(0, 1) random variables. Let us assume that N > n, 
and N — n = v is fixed. Consider a nonnegative-definite random matrix 

M n , N = AA l . (24) 

The ensemble of random matrices M n ^ is known as the real Wishart distribution in the statisti- 
cal literature or the Laguerre ensemble in the mathematical physics. The empirical distribution 
function of the eigenvalues of M n ^ converges to the Marchenko-Pastur law as n — > oo ([8], [IT]). 
The density of the Marchenko-Pastur law is given by 



^V^x : if < x < 4, 



pmp(x) = \ ; • - ' (25) 

| (J : otherwise. 

Our goal is to study the distribution of the smallest eigenvalues of M n ^ in the limit 
n — > oo, N — n = v. It can be shown (see e.g. [4], [5], [3]) that the smallest eigenvalue of M n ^ 
are of the order of n~ 2 . Moreover, if we consider the rescaling at the hard edge of the spectrum 

Jn,N) 

Af^ = L x< z < n , (26) 
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one can show that the random point configuration {Q , 1 < i < n} converges in distribution 
on the cylinder sets to the pfaffian random point process on (0, +00). The /c-point correlation 
functions of the limiting process are of the same form as in ^j, where K(x,y) is again a 
2x2 matrix- valued kernel. The formulas for the entries of K(x, y) are similar to ([5]) with the 
important difference being that the Airy kernel KAiry(x,y) is replaced by the Bessel kernel 

R , v _ Ju(>/x)y/yJv+\(y/y) ~ Jy{^)\fxJ v+ l{^X~) 

Bessel(X,y) 2(x - y) ' 

where J v (x) is the usual Bessel function of index v. 

Define ^ 

G n , N (t) = ^M n ,jv + ^j , (28) 

where t is a real number. Then 

9k(t) = 9n,N,k(t) = n~ 2k TrG k n<N (t) = n~ 2k Tr (m^n + J) = f^'^ + t 2 )~ k (29) 

is a "local" statistics for any positive integer k = 1,2,3,.... Indeed, one can show that the 
eigenvalues from the bulk of the spectrum give vanishing contribution to gk(t). In particular, if 
£ = te> 2 E N} is a random point configuration of the limiting pfaffian process, then 

r+00 1 

+ t2 y k = I j^rw piix)dx < °°' (30) 

We are interested to study the joint moments of the linear statistics gk(t), k > 1. Let, as 
before, K = (ki, . . . , kj), j > 1, denote a multi-index, and stand for the corresponding joint 
moment 

3 

m K = EH( gi (t)) k <. (31) 
1=1 

The following theorem holds. 

Theorem 1.3. Let K be a non-zero multi-index, then the following equation holds: 

1 v— ->. 2 v— -\ 

(v - t~ 2 + -)m K+ei + m K +e 2 + m K+2ei + 2 ^ lkimK- ei +e l+2 - 75 I m Jf-e,+e I+1 

i=l 1=1 

11 

= -T^niK- (32) 
A/so, £/ie following "boundary" condition holds: 

(v-t~ 2 -\ )m e , + m e , + m 2e , = t~ 2 . (33) 

n 
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We recall that v = N — n > is the difference between the dimensions of the rectangular 
matrix A n ^. 

We finish the Introduction by the discussion of the complex Wishart ensemble. Let A = 
A n! N = ~^( a ij) be a rectangular n x N matrix with {SRa^-, Qa^, 1 < i < n, 1 < j < N,} i.i.d. 
iV(0, 1/2) random variables. As before, we us assume that N > n, and N — n = v is fixed. 
Consider now a nonnegative-definite random matrix 

M n<N = AA*. (34) 

The ensemble of random matrices M n ^N is known as the complex Wishart/Laguerre ensemble 
of random matrices. Consider the rescaling of the eigenvalues at the hard edge of the spectrum 

Jn,N) 

Af^ = L, x<,< n , (35) 

It follows from the results of [1] that the random point configuration {£■ , i > 1} converges 
in distribution on the cylinder sets to the determinantal random point process on (0, +oo) with 
the correlation kernel given by the Bessel kernel (|27p in the limit n^oo, N = n + v. Define 
Gn,N, 9k(t), k>l, and txik in the same way as in (|28p . (|29|) . and (|31|) . The following theorem 
holds. 

Theorem 1.4. Let K be a non-zero multi-index, then the following equation holds: 

1 j 1 j 

(y + -)m K + ei + m K+2ei + ^ lkim K - ei+ei+2 - ^ lkim K _ ei+ei+1 

n i=i i=i 

11 

= --nm K - (36) 
n t z 

Also, the following "boundary" condition holds: 

iy + -)m ei + m 2ei = t~ 2 . (37) 
n 

We recall that v = N — ra > is the difference between the dimensions of the rectangular 
matrix A n jv- Theorems 1.3 and 1.4 will be proved in Section 3. 



2 Proof of Theorems 1.1 and 1.2 

Let us start with the proof of Theorem 1.1. Our first goal is to establish (I17p . To this end, we 
consider n 1//3 m ei = E (ra -1 / 3 (n + TrG n {z))} , where, as before, G n (z) = (A n — 2 — zn -2 / 3 ) -1 . 
By using the resolvent identity 

G n (z) = (A n -2- zn~ 2/z Y l = -(2 + zn~ 2, *)- l Id + (2 + zn' 2 ^)" 1 A n {A n - 2 - zn" 2 ^)- 1 = 
-(2 + zn- 2 l z Y x ld + (2 + zn- 2 ^Y x A n G n , (38) 
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we arrive at 



n 1//3 m ei = n 



,2/3 _ (2 + zn -2/ 3) -l n 2/3 + (2 + zn -2/3 ) -l n -l/3 E £ AijGji. (39) 



Here = ^% denote the matrix entries of A n , and Gij denote the matrix entries of G n {z). 
To calculate EAijGji, we recall that random variables Aij, 1 < i < j < n, are independent. 
Therefore, we can first fix all matrix entries (up from the diagonal) except A^ and integrate 
with respect to A^. Applying the Gaussian decoupling formula 

E V f( V ) = a 2 Ef'( V ), V ~ N (0,a 2 ), (40) 
with rj = A^ and f(rj) = Gij, and taking into account that Var(Aij) = 1+ ^ J , and 



9Gh = f -GkiGji - GkjGu : i / j 
dAij \ -GkiGji : i = j, 



(41) 



we arrive at 



,V3 



n 1 m ei = n 



i 2 l 3 - (2 + zn" 2 / 3 )- V/ 3 - (2 + zn- 2 ^)- l n-^E ^ {GjiGji + GuGjj) ■ (42) 



The term n 4 ^ 3 E^^- GjiGji is equal to 

" 4 / 3 E J2 GjiGji = En~ 4 / 3 Tr(G 2 n (z)) = m e2 . (43) 



n 



To deal with the term n" 4 /^^^ GuGjj, we rewrite it as 

n" 4 / 3 E GuGjj = n" 4 / 3 E ((TrG„(z)) 2 ) = n~ 4 / 3 E ((-n + n + TrG n (z)) 2 ) (44) 

= n 2/3 - 2n~ 1/3 E(n + TrG n {z)) + n~ 4/3 E ((n + TrG n {z)) 2 ) = n 2/3 - 2n 1/3 m ei + m 2ei ■ 
As a result, we obtain 



n l ' 3 m ei = n 2 ' 3 - (2 + zn 2 l 3 y x n 2 l 3 - (2 + zn 2 l 3 y x m^ - (2 + zn 2 ' 3 )^ 



n 



2/3 



+2(2 + zn' 2,3 )' 1 n 1/3 m ei - (2 + zn- 2/3 )~ 1 m 2ei , (45) 
which is equivalent to 

n 2/3 (1 _ + ^n-2/3/2)- 1 ) _ m2ei ( 2 + zn" 2 / 3 )" 1 - m e2 (2 + zrT 2 l 3 y l 

^(^(l + zn- 2 / 3 ^)- 1 ). (46) 
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After trivial arithmetical simplifications, this leads to (|17p . The formula (|16p can be proven 
along the same lines if one starts with n'^m^j. One can say that the formula (|17p gives 
us the boundary term in the recursive system of linear equations satisfied by {mx} since it 
corresponds to K = 0. Turning our attention to (|16p . we write 



n 



E 



1/3 



m K+ei = n 1 /^ « 1 (z)) fcl+1 X\(9n,k^)) h 



1=2 



n 



1/3 ( n + TrG „) L- 2 /3( n + TrG n )) fcl Yl (n- 2l / 3 TrG l 



l>2 



we then rewrite, as before, the first term n 1 < 3 (n + TrG n ) as 

^ 3 (n + TrG n ) = n 2/3 - (2 + zn" 2 / 3 )- V/ 3 + (2 + zn" 2 / 3 )- 1 /^/ 3 £ A^-G 



?i 



(47) 



(48) 



This leads to 

n 1//3 TOK +ei = n 2/3 m K - (2 + zre" 2 / 3 )™ 1 ?! 2 / 3 ?™^ + 



(2 + zn- 2 / 3 )- 1 n- 1 / 3 E 



(£ AijGji) (n- 2 / 3 (n + TrG n )) * J] (V^TrG* 



2>2 



(49) 



As in the case K = considered above, we fix all matrix entries (up from the diagonal) except 
Aij, and apply @Q) with ry = and /(r?) = (n~ 2 / 3 (n + TrG n )) kl ]Ji>2 {n- 2l/3 TrG l ) k ' . 
Taking into account (|4ip and the equation 



OA 



-2l(G l 



(50) 



one then obtains (|16p after some simple algebraic calculations. Theorem 1.1 is proven. 

The proof of Theorem 1.2 is quite similar. The only alteration required in the GUE case is 
that one needs to replace (f4*T]) with 



8G. 



and 



dRe(Aij] 
dG kl 



-Gki.Gji — GkjGu : i ^ j 
—GkiGji : i = j, 



-i(G k iGji - GkjGu) for i ^ j. 



dIm(Aij) 

The remaining calculations are very similar and are left to the reader. 



(51) 
(52) 
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3 Proof of Theorems 1.3 and 1.4 



The proofs will be similar to the ones given in the previous section. Let us start with the 
proof of Theorem 1.3. Our first goal is to establish (|33p . To this end, we consider n _1 x 
m ei = E (n~ 3 TrG nj N(t)) , where G n ^{t) was defined in (|28j) . By using the resolvent identity 



G 



we arrive at 



n,N 



n m, 



{M n>N + 



n- 



1 1 



n 

l 2 



n 



Id - ^AA f G 



n,N 



ei 



t 2 



n 



t 2 



E A ip AjpGji. 



(53) 
(54) 



l<i,j<nl<p<N 

Here Ai p = 2% denote the matrix entries of A n , and Gji denote the matrix entries of G n ^{z). 
To calculate EAi p Aj p Gji, we again use the Gaussian decoupling formula (HOI) and the equation 

dG kl 



-Gki{A t G) p i — (GA) kp Gi 



(55) 



Therefore, 



n l m ei =- + —-Ej2 E A iv ( G A GA )i P + Gu(GA) jp ) - — -E ^ £ G « 

l<*,i<™ 1<P<^ l<i<nl<p<AT 

4 + 4t^E [Tr((L44'G) + Tr(GAA*)TrGl - ^^ETrG. (56) 



t 2 

Using the identity GAA 1 = Id — ^G, the last formula can rewritten as 



n Trip 



n + v 

t 2 me l' 



(57) 



which implies (|33|) . The formula (|32p can be proven along the same lines if one starts with 
n~ 1 mK +ei . Let us write 



n rriK+a = n 



~'E (<&(«)) flfoC*))' 



l=i 



E 



n- 3 TrG n , N )H(n- 2l TrG l r 



T n,N 



k, 



1>1 



(58) 



Using the resolvent identity, we can rewrite the first term in the product as 

Tip 1 TlP 

G n ,N = ~ ~p AA G n ,N- 

After integration by parts and a few lines of careful calculations, we obtain (|32p . Theorem 1.4 
can be proven along similar lines. 
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4 Non-Gaussian Case. 



The generalization of the Gaussian decoupling formula (|40p to the non-Gaussian case can be 
found, for example, in [7J: 



d k f, 



dx k ® 



+ e, (59) 



k=0 

where £ is a real random variable such that E (|£| p+2 ) < oo, q, Z > 1, are the cumulants of 
the random variable £, complex-valued function f(x) has first p + 1 continuous and bounded 
derivatives, and the error term satisfies the upper bound |e| < B] t sup x \^-^(x)\E (|£| p+2 ) with 
the constant depending only on k. 

It is conjectured that the distribution of the largest eigenvalues in Wigner random matrices 
is universal provided the fourth moment of the matrix entries is finite. Currently, we are unable 
to prove this conjecture. Instead, we speculate below on the possible approach to extend the 
results of Theorems 1.1-1.4 to the non-Gaussian case. Let us consider a real Wigner random 
matrix A n = -^(oij)™j =1 , and assume that the entries (a^ = ajj)i<j are i.i.d. centralized 
random variables with the unit variance and the finite fourth moment. In addition, we assume 
that the diagonal entries an, 1 < i < n, are i.i.d. centralized random variables, independent 
from the non-diagonal entries. We assume that the diagonal entries also have the finite fourth 
moment. Let us also assume for simplicity that Var{an) = 2. In an attempt to extend the result 
of Theorem 1.1 to the non-Gaussian situation, we apply the generalized decoupling formula (|59p . 
To be specific, let us concentrate our attention on the "boundary" equation (fT71) . Looking at 
HMD, we apply d59j) to EAijGji. Since ci(ay) = 0, c 2 (aij) = 1 + (%, 03(0^) = c 3 , for i < j, 
and Eafj- < 00, one might wish to truncate (j59[) after the first three terms (i.e. p = 2) to obtain 

n -i/3 E AjGji = -n- 4 %^(^ + GuGjj) + n - n / 6 E ^(SG^GaGjj + G%) + e. (60) 

ij ij ij 

The first sum in (|60p is the same as in the Gaussian case. The hope is to show that the second 
sum and the remainder term give negligible contributions in the limit n — * 00. However, it is 
currently unclear to us how to efficiently bound the terms n~ 11 / 6 E^ ij - GijGuGjj and 

n -ll/6 x 
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